





(tt: 


HEN 
Di 


1 
1 





Ley 
TM 








SC Vë 


32: t 
EE 


1 eae 














سا 


—— 
— 
we 
re 
L 
u 
D 
—— 
DH 
س‎ 
س‎ 
— 





Lupa Thco c c E M. uM L| 


qun Fir : 
BER: 


The Engineering Control Problem 


Control problem can be stated as follows: 
A physical system or process is to be accurately controlled in order that certain control criteria, or specifications, be 
satisfied. | 
In general, solving a control problem involves the following steps. 
° Aset of performance specifications is established ( requirements). 


e As a results of the performance specifications a control problem may be exists. 
To solve the problem and additional steps are required ; 


° Aset of differential equations ( Modeling) that describe the physical system is formulated. 
° Using the classical control-theory approach to analysis the system results (Simulation). 
° Accordingly, Control system Design must be implemented to improve the system 


° Using the modern control —theory approach , the designer specifies an optimal performance index for the 
system ( 1.6 the designer yields the necessary structure to minimize the specified performance index, thus 
producing an optimal system. 


° [f the physical tests are unsatisfactory, iterating the above steps. Problem formulation 
e a è 5 S Mathematical 
o The process can be summarized as in the following schematic diagram modal of 
Physical Conceptual 
system aspocts 
solution of 
mathematical 
Solution translation nw 


System modelling 


2 M od e | | | ng: is the description of the behavior of physical system by mathematical approach. 
° These mathematical equations must be 
— Not too simplified 


— Not too complicated. 
° The results of these models must be reasonably accurate. 


° Types of Models 
e  Asystem model can be classified into 


— Static Model : 


(1) the system whose some of its elements behavior fast relative to other elements is often 
modeled as static model. For example , the switch time of a thermostat is fast compared to 
the time required for the room temperature to change appreciably. Thus, the thermostat is 
modeled as a static model. 

(2) the system whose output at any time depends on only the input at that time. For example, 
the resistance may be considered as static element, because its current depends on only the 
voltage applied at present and not on past voltage. 


— Dynamic Model : 


it is the model whose present output depends on past inputs. For example , the present 
position of a car depends on where it started and what its velocity has been from the start. 











System model classifications 


system models are classified according to the types of equations used to described them. 
The following family tree illustrates the major system classifications. 


Classes of | 
systems 
Distributed geen 
parameter parameter. 
Stochastic Deterministic 1 
Continuous Discrete 
: time time 
cJ 
Time Constant 
` varying coefficient 


Major classes of system equations. 
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Modeling steps 


An outline of the analytical approach to model physica systems can be represented by the 
following steps : 


^ Purpose of the 
model 
p Define boundaries i 


Gef Postulate a structure 
à which mimics observed 


or desired effects ` 





Pd Specify variables 
of interest , 


Develop mathematical 
description of each 


modei element 





Apply physical laws 
of continuity and 
compatibility 


Manipulate 
equations 


Modify model 
if necessary 


f 
| 


ام موس 


Final form of 
mathematical model 


Analyze. compare 
with real world 
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(b) Steps in modeling 


Linear systems (summary) 


1. In practical most physical systems are non-linear. They are often approximately modeled by 
linear equations mainly for mathematical simplicity. This simplification is Satisfactory as long as 
the resulting solutions are in agreement with experimental results. 

2. For linear systems the concept of a transfer function can be used. 

3. Alinear systems in state space can be represented by the form x = Ax ~ 

4. Linear systems are a systems to which the principle of superposition applies. 

5. A linear systems can be analyzed in time or frequency domain. 

6. For linear systems a sinusoidal input leads to a sinusoidal output with the same frequency. 


7. There is only one equilibrium point in linear system xx = (0 | 


8. The stability concept depends on the system poles locations . 





Non-linear systems 


* The non-linear systems differ from linear systems greatly in that 
the principle of superposition can not be applied. 
* . Non-linear systems have many phenomena that can not be seen in linear systems. 


such as : 


1. the dependence of the non-linear systems behavior upon the magnitude and type of the 
input. (LE the nonlinear system behave completely differently in response to step inputs of 


different magnitudes. 
2. Stability of the system is no longer just a simple function of system poles location 
but depends on the number and position of system Equilibrium points. | 
3. The correlation between transfer function pole and zero locations and time 
response behavior is generally invalid. 
4 Finite escape time: The state of an unstable linear system goes to infinity as time 
approaches infinity ; a nonlinear system's state , however, can go to infinity in finite 


time 


For many nonlinear systems which have sufficiently small and smooth nonlinearities can be treated 
by using equivalent linearization techniques such as a linear approximate model, obtained by 
linearizing about a known nominal operating point. 
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Wight control ( open loop) 


Wight control ( close loop) 


An aircraft's attitude control ( roll angle stabilize) 
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° Control of the speed of a turntable 
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Problems 


An automobile driver uses a control system to control steering ofthe car on the desired 
direction Sketch the block diagram for this system. 





A control system shown used a human operator as part of a closed loop system . Sketch 
the block diagram for this system. 





RN 











3. The student —teacher learning process is inherently a feedback process intended to reduce 
the system error to a minimum. Construct a feedback model . 


4. Draw a block diagram representation of a thermostatically controlled electric oven in the 
kitchen of a home. 


5, An elevator — position control system is used in an apartment building . Draw the block 
diagram representation of the elevator in a three-floor building which obtains the desired 
floor reference position as a voltage from the elevator passenger pressing a button on the 
elevator, and compares this voltage with a voltage from a position sensor that represents 
the actual floor position the elevator is at. The difference in an error voltage which is 
amplified and connected to an electric motor that positions the elevator car to the 
desired floor selected. 


Solution 


y. Desired Actual 
3 


floor floor 
position ` Elevator 
dynamics 
Floor position 
Sensor 


OV, OF Ok, 
Floor 1 2 3 


















6. The automatic depth control of a submarine given below is an interesting contro! problem. 


The desired depth is set as voltage from a calibrated potentiometer. The actual depth is 
measured by a pressure transducer. Draw the block diagram representation for this system 
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6. For antiaircraft radar tracking control system shown below . Sketch the block diagram for this 
system e 
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Mathematical models of 
physical systems 


Most of physical systems are characterized by the term mathematical model 


Mathematical model is defined as a set of equations that represents the dynamics 


of physical systems. 
To find mathematical model the following informations must be provided 
1. Appropriate basic physical laws. 
2. The system elements. 
3. The way the system elements are interconnected 


Linearizing the non-linear systems ( for certain case only). 


Electrical systems 


e The basic system elements are 


Elements model Model equations 
TUIS EM /-I* 
٣ کا کا‎ $^ 


| e Resistance (R ) 


| © Inductance (L) 





e Capacitance ( C ) بب‎ es Ke 


° The basic physical laws 
e Ohm's law. 
e Kirchhoff node law 


e Kirchhoff loop law 
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Examples 


For the following systems find the mathematical model 
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3. Forthe following inverting amplifier find the mathematical model. 


The equation for this circuit can be obtained as follows: Define 


B= . do 


h = سس‎ -— 


l 
’ 2 





Since only a negligible current flows into the amplifier, the current ij must be equal to 
current i, Thus : 











CC 
Since K(0 — e’) = e and K> 1, e’ must be almost zero, or e’ = 0. Hence we have تت‎ 
~ 1 
d Wi 7€, | 
M Om | | 
01 | E 
, | 
d 
g p. i 
0 R; l i 


Thus the circuit shown is an inverting amplifier. If R; = Ry, then the op-amp circuit 
shown acts as a sign inverter. 
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Mechanical Translational 
system 


e The mechanical T.S elements involved are 
e Mass | TET f(t) | DU 
* Linear spring 
e Damping (friction ) 
e The basic laws are 
Newton's law 


e  Hooke's law 





EE M (f) 


fit) | f(t) ft) 
. Mass "T Friction Spring 


f(t) = Ma(t) =u = اوق‎ f(t) = ې سوا چا‎ II - ٨-٧٧ 


Examples 


e For the following simple systems . Write the mathematical model. 
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For massa 


Mechanical Rotational system 


The mechanical Rotational system elements involved are 


* Moment of inertia E "UE Bal : Hole 
Where t(t)( T(t)) is applied torque, f=) P dr d 


J is moment of inertia 


O is angle of rotation 


* Torsion spring prz IE 
K is the spring coefficient | Tli) z K BUI = LI 
t $4 95 4 د‎ Y i jer Fe 
° Damping (friction ) Iz T ] 19 Hi) E 
B is the damping coefficient 33H 55" Tf} = EH = b GU) - th (f)]| 
O4 T 95 ili d / 


The basic laws are 
° Newton’s law 


Hooke’s law‏ ه 





Conversion between Translational and 
Rotational Motions 


e In motion-control systems it is often necessary to convert rotational motion into translation . For 
example a load may be controlled to move along a straight line through a rotary motor-and 


screw assembly Or through a rack-and-pinion Or through a pulley by rotary motor . The systems 
given below can be represented by simple equivalent inertia 
e Where Lis defined as the linear distance that the mass travels per revolution of the screw 
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Examples 


For the following systems find the mathematical models 





T(t) يرم‎ 


(a) 


+ ې سو‎ KkO(t) 


dO(t) 
1) - J— +B 


POA) 
dt? 


Electro-mechanical systems 


The most common electro-mechanical system are 


Air gap 





— Dc motor Ue 
Motion شا‎ Kc) 
Dc motor armature control ( constant field current) | = : | E 
á - Corc 
E,-L,— 
dt 


* 


E 2 Ri 4L, Se X p, 
01 


40 

E, = k, — 
b b dt 

Flux p=k, *i, 

Torque بپ‎ a 0 j L = لک‎ 

Apply Newton's second law of motion yields 
2 

d'8 LB dé 


T=J——+B— 
dt dt 


Dc motor field control ( constant armature current) 
° Drive the mathematical model for this type. 





potentiometer 


A potentiometer is an electromechanical transducer that converts mechanical energy into electrical energy 
The input to the devices is in the form of a mechanical displacement. 


When a voltage is applied across the fixed terminals of the p: entiometé&) the output voltage , which is 
measured across the variable terminal and ground, is proportional tokinput&isplacement. 


e(t) = k,0(f) 
where E 
0 


bo, 


N = no of turn 
227٧۷ 





e(t) 


For more flexible arrangement is using two potentiometers connected in parallel 





e(t) = K,[8, (4) - 8,0) 


e(t) 


Tachometers 


e Tachometers is device to convert mechanical energy ( speed) into electrical 
energy. 


e The dynamics of the tachometer can be represented by the equation 





Gear system 


The gear system can be modeled under the following conditions 
1. The gears are rigid. 
2. Gear teeth are meshing perfectly. 
F, Gears have no inertia. 


Then the distance traveled by one gear must be equal to that traveled by other, or 


The force exerted by one gear must equal the reaction force of the other at the point of 
contact : 





Antenna system 


e For the following Antenna system find the schematic diagram 
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Lever mechanism 


e Lever system is shown as 





e For small angles from horizontal , the total motion y equal the sum of motion due to x1 with x2 = 
0 and that due to x2 with x1= 0. 
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Liquid Level system‏ ھ 


e Consider the system given below. In steady-state 


net volumetric in flow rate= change per unit time of the volume | 


y 


Gm Cha E 


í 
where v is tank volume 





The relation between H and Qo is 


where R is Fluid resistance 
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e  Fortwo-tank liquid system 


Tank | Tank 2 
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Modeling an Inverted Pendulum 


The schematic diagram of inverted Pendulum is given as 





let's assume that 

M mass of the cart = 0.5 kg 

m mass of the pendulum = 0.2 kg 

b friction of the cart = 0.1 N/m/sec 

L length to pendulum center of mass = 0.3 m 
| inertia of the pendulum = 0.006 kg*m^2 

F force applied to the cart 

x cart position coordinate 

theta pendulum angle from vertical 


Force analysis and system equations 


e Below are the two Free Body Diagrams of the system. 





bx 
Hiction 





[ -x 


— — Be xX 


e Summing the forces in the Free Body Diagram of the cart in the horizontal direction, the 
following equation of motion is obtained: 


Mx+bx+N=F 


e Note that you could also sum the forces in the vertical direction, but no useful information 
would be gained. 


e Summing the forces in the Free Body Diagram of the pendulum in the horizontal direction, 
an equation for N is: 


N=mi+mldcosé -mlé“siné 


These set of equations should be linearized about theta = Pi. 


e Assume that theta = Pi + @ ( ọ represents a small angle from the vertical upward direction). 
Therefore, 


cos(theta) = -1, sin(theta) = -@ , and (d(theta)/dt)^2 = 0. 
e After linearization the two equations of motion become (where u represents the input): 


| I+ ml* lé -melo = mix 


M 4 më +bx- mld =U 


Model a belt drive printer 


e A schematic diagram of a belt drive printer is shown as 


Ti 
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Printing 
device 













v(t) - 
Printing Controller Sensor 
: a avi : : 
1 = eO — U) = KAN 
device P w= اس‎ dj (LS Ku y 
position 


e inthis model , a light sensor is used to measure the position of the printing device. 
e The belt tension adjusts the spring flexibility of the belt. | 
e Where 


dv, 
^ dt 


y, —- Ay, Qu 


EE x 


System model 


The relationship between the pully rotation ( Op) and printing device displacement y is 








y=r0, 

The tension T1 is L =k(r0-ré, (- ٧) ۳6۷ — y) 
The tension T2 is بر‎ E k(y E re) 

d^y - 
Then the total net tension is Í; 0 I; =m PE = ZKk(r8 - y) 
In motor side ( assume L=0 ) then 
The field current is Ee 

رر 

Motor torque Tm is T, = k i, =f 22 


Hî 
R 
The motor torque provides the torque to drive the belts plus the disturbance or undesired load 
torque, so that Tm=T+Td 
Where 


d'O d 
T=J +b 4H(T -T 
d? di ec 
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Thermal system 


Thermal system consists from 


Thermal resistance ( Rt) 


A wall of area A separates regions with temperature T1 and T2 


The heat flow rate q, in units of heat per unit of time, is proportional to the temperature 
difference T1-T2 and to the area A and flows toward the lowest temperature . 


The constant of proportionality is the heat transfer coefficient h . 


A 
q o FA —T,) 





kA L-T 
(B ~T,) = hA(T, E 


gegen 
tegen 
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Thermal Capacitance 


let q be the net heat flow rate into a volume V of a material with mass 
density p and specific heat c ( = heat required to raise the temperature of a 
unit mass by 1°). q 


p 


MP 


The inflow q of heat per second must equal the change per second of heat 
stored in volume V. 


q oke = CT 


Where C, is called a thermal capacitance. 


Examples 


o For the following thermal system find the mathematical model. 





° The behavior of the system can be modeled as 
dT 
E The net flow is q;—q,-C,.— 
dt 
Z 
2. The output heat is q, = — 
R, 
3. Therefore the mathematical model is 
dT 
Ja =1 ې‎ R,C, nw 


dt 


Model a liquid-in-glass thermometer 


To obtain a linear model relating the thermometer indicated temperature, T2, to the measured 
temperature, T1 


According to Fourier's law, the heat flow rate through the glass bulb is.proportional to the 
temperature gradient within the bulb, | 
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q = kA 
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area,A wä ) 


Mercury, mass rr 


ef 88s See 
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temperature, 74 Bulb thickness, تر‎ 


اس 
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Where k is thermal conductivity of glass bulb and A its surface area. 


The temperature difference across the glass bulb will be (T1-T2), B bulb thickness 


dT T1—T 2 
dx B 


ri-f2 
B 


Net heat flow = Rate of change of internal heat C = 





q = kA 
dH 
41 


where H is the internal heat of the mass m,of mercury at uniform temperature, 
T2 is given by 
H =m cp T2 


where cp is specific heat. dT 
- kAT , + mCpB m = KAT. 
dt 


49 


A temperature-control system 


> This example is an industrial 5 control . Draw the block diagram 
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Fluid system 


Fluid system elements are defined as 
Control valves 
Two common types of valves are shown below 
— The supply pressure is p.. 


— The output pressure is p. 
— For x=0 the output port is just blocked off. 
e Valve flow increases with both x and valve pressure drop ( ps- p). 


q=k,x +k (P,= P)=KxX+K, Pi 
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2- Pneumatic tank 





° The flow q through R ( Fluid resistance) can be written as 
C 
. Due to this flow the tank pressure will raise according to 
d 
= C P o 
dt 
۰ Rearrangement the above equations yields 


p 
= — . 
RC 7 * د‎ = Pi 


: So that the transfer function is 


p, l 
p, RCs+l 





3- Pneumatically actuated valve 


Pneumatically actuated valve is used extensively in level control systems. 


Diaphragm 
area 


The relationship between control pressure pi to pressure p is represented by 


dp 
eC. 
17 


The block diagram is given as 


Di 





Actuator 


P:i P 


q = R 


Due to this flow the pressure will raise according to 


Diaphragm motion x is to small that the capacitance C of the space above it is about constant 


C i 
LA 


The downward pressure force (Ap ) on the diaphragm must be counterbalanced by the spring 
force ( kx ) | 


Ap = kx 


the model for the valve flow q, is given by dq. 3 kx T k pd 
Where pd is the pressure drop across the valve. 


Finally ,the following transfer functions can be derived 


wA 
ps) — 1 Qs) ^ 7k 
p(s) RCscl p(s) RCs«l 
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If a load consists of mass ( m ), spring (k) and damper ( b ) is connected to the hydraulic cylinder as shown 
below. Then 


the pressure on piston is p, and the force balance equation is 
p4 = my + by + ky 
The valve flow is Oz k x سل‎ k, (D, ES p) 


The flow q to the cylinder with compressibility flow associated with pressure variations 5 


F, 

q = Ay + — p 
D 

V 


Where V is the volume under pressure, and compressibility flow is — p 
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Linearization of non-linear system 


The linearization procedure is based on the expansion of the nonlinear function into a Taylor 
series about the operating point and this can be done by assume that the variable deviate only 


slightly from some operating condition i.e neglect higher-order terms of Taylor series. 


Steps for linearization 


Define the nonlinear element. 
Write the nonlinear equation which represents the nonlinear element. 
Assume that the variable deviate only slightly from some operating condition. 


Approximate the the nonlinear element by using Taylor series . 


Let the input — output model for a static element be written as 


y = f(x) 
If the normal operating condition corresponds to xo, yo then 


eeooed 


ach 2 3 
f(x)= f+ Gx) oF 
on ١ & 





X=X, 
if the variation x — xo is small , then the above equation can be written as 


y =f(x) =f(xo) +K (x-xo) 


k 2 S. 
ax 


Where 





2-0 f(x.) — — — A 
Equilibrium 
Operating point 


Then the linearized equation can be written as : 


y -yo-k(x- xo) 


Example 1 find the linearized model for the following system about the operating point 


f(x)=Scosx x, = Z 


-—3 sin(x) _ = =—! 
ER: " 
0 


F(x,) Seet 
y-X0)2 k(x —x,) 


f(x)- 02 -5(x- 77) 
f(x) = e 74) 

















E le2 : _ 
INT 8 د‎ RET, r= 74 
at” dt 
e Nonlinear portion f(x) = COS x 
af | 42 | 
k zl az = سب‎ 5111 X — ——— f (xq) =cosx|._ 7 
X odi 4 
Ep -x-7, ét x= +4 
f (x)= X xa 
dE: 7/ 
Linear model er) deem) 4 کل‎ 
e 
dt dt 2^ 2 
016) dër 42 2 
017 ^ o 


E کا یی‎ 
dt 2 


Example 3 
V, (s) 


Find the transfer function dor the electrical network shown in figure below, which contains a 
nonlinear resistor whose voltage-current relationship is defined by 


° : 0.17, 
l= يره‎ 


Solution 
Use kirchhoff's voltage law to sum the voltage in the loop to obtain the nonlinear differential equation. 
The voltage across the non-linear resistor is 1 
The current-voltage equation around the loop is د‎ 210 AR 

L s. -101In ا‎ 20 = v(t) oOo o AO 
Evaluate the equilibriu solution, by set the source 
Voltage signal equal zero i.e 

with v(t) = 0, the circuit consists of a 20 V battery in 

Series with the inductor and non-linear resistor. 
In the steady-state, the voltage across the inductor will 
be zero, and a given battery voltage will be cross the 
resistor, then 





Nonlinear 


d' | 
resistor 


aim, ¥ 





= 3* 4 EE, a aru 
Thigcurréat is RE equilibriátft value, aad the total current is 


tø‏ د 


Then the non-linear part is =) — E 
Linearize this part yields ` | - 


دو الا 
se‏ اط اکا 


di 1 


o و هس‎ lg 
^; ` dä 


The linearized equation then is L —— +10 In) T Ze — 20-2 v(t) 
dt 2 i 


0 





Letting L = 1 and io = 14.78 yields + 0.67761 = v(t) 
1 
| V (s) 
Taking the Laplace transform with zero initial conditions yields Os) m ———— — 
s + 0.677 
d A dói 
The voltage across the inductor about the equilibrium point is Vi (t) = b Me + ói) = Ne 
Taking the Laplace transform, 
. . sV(s 
V, (s) = Lsói(s) = sói(s) = 21. 
ووا‎ 0.677 


Then, the final transfer function is 


V, (s) D S 
Vis) s+0.677 


(x, —x,(0)) 


i- 
e 

f 
bd 


2 





ix E 


p= mRT 


rar, و‎ d 1 


oP 
OV T= 


(x, -x(0)) +L 


Multi-variable 
The function is CHEN 


Taylor series expansion is 


fex.) = A (0}x4(0) +2 


1 





X =x, (0) 


Example consider the perfect gas ۷ 














OP oP 
Ps P+ (T-T,)+— 
OT 7 =1) OV ۶ 
H zu dis df 
GEI _ mR = mR 
OT \T=% PF اا‎ F, 
y -V, y -V, 





op = P-R; =T- hð, =V -Vo 





x=0 X-—J 
Where 31101-02٥۳6 5 
2 
X 5 ES 60S, - 


X, = X, + (x, T L)x, +, BIR E, 


¥+4x° +(x? +I} =0 x=0 


Problems 


2 
E PR LUN, C 


1) 
| ar dt 





2 

(2) Q = KA ^E (p, — Pa) 
P 

(3) Linearized the following system 


(5) Linearized the following system 


(6) Find the linearized transfer function for the electrical network shown below, if WEEN 
vit) 
ځ د و‎ A) Nonlinear =, LE 






resistor 


$ 
| 
f 
d 
i 





(7) Find the linearized transfer function for the electrical network shown below, if 





Nonlinear | pa 4.2, 
a | vtt) zu 2l LU) 
resistor | ` 
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1- Block diagram 


Block diagram: is a representation of the function performed by each component and of 
the flow of signals. Such a diagram depict the interrelationships which exist between the 
various components. 


In a block diagram, all system variables are linked to each other through functional blocks. 


Functional block is a symbol for the mathematical operation on the input signal to the block 
which produces the output. 


Transfer function: (T.F) 


— This phrase is very often used to characterize the input-output relationships of linear 
time-invariant systems and it is defined to be the ratio of the Laplace transform of the 
output to Laplace transform of the input under the assumption that all initial conditions 
are zero. 


Consider the linear time-invariant system defined by the following differential 
equation: 
(m) (8-1 : 
agy + وړچ‎ occ db duy d du 
m (n—1) 1 . 
= Bx + b,x + =e +h ct bx (n > m) 
where y is the output of the system and x is the input. The transfer function of this 
system is obtained by taking the Laplace transforms of both sides of Eq. d 
under the assumption that all initial conditions are zero or 
: : YD bas” + Bsmt p °< ګر سل‎ + b, 
Transfer function = G(s) = nn = p Tas on ټپ‎ h د‎ ya cca, gs TTT 


The transfer function of a system is a mathematical model in that it is an operational method of 
expressing the differential equation that relates the output variable to the input variable. 

The transfer function is a property of a system itself, independent of the magnitude and nature 
of the input or deriving function. 


T.F contains information concerning dynamic behavior, but it does not provide any information 
concerning the physical structure of the system. 


If T.F of the system is known , the output or response can be studied for various forms of input 
with a view toward understanding the nature of the system. 


If the T.F of the system is unknown, it may be established experimentally by introducing known 
inputs and studying the output of the system. 











Transfer 
function 
G(s) 


Ys) 


Closed-loop system subjected to a disturbance 


When two inputs ( the reference input and disturbance) are presented in linear closed-loop 
system , each input can be treated independently of the other. 


In examining the effect of the disturbance N(s) , the input reference R(s) is assumed zero then 
the response due to disturbance Cy (s) can be found from 


On the other hand , the response due to reference input R(s) can be found assumed zero 
disturbance input such as 


The total response C(s) then is obtained by adding the two individual responses. 


C) Gu — 
Ns) | t G,(9)G,(9H() 





1 Disturbance 


Cs) ` G,(5)G,(s) 
Bai 1+ G,(8)G,(s)H(5) 
(5) = Cs) +04 





D 
, 


"ان 
lee‏ 


/ کې‎ 
ay SE 


Example 


For the level tank control find the over all transfer function. 


Tank | Tank 2 "ët: 





eege mc 


Problems 


Find the over all transfer function of the Dc motor armature control . 
For the following system show that the over all T.F is given by 


E E 
lC RCs + (RC) + Roly + RO 





en‏ ند ma orti‏ حون 3 NIT‏ لد Jg)‏ مره ich‏ ېښ RR E‏ ولوار 


2-State-space 6 


State is the smallest set of variable( called state variable). Or is defined as minimum 


information required at time to, to find the response to subsequent inputs( t > to) 


State vector is a vector that determines the system state X(t) for any time t 2 to. 


State space. The n-dimensional space whose coordinate axes consist of the x1 axis, x2 
aXIs,............, XN axis is called a state space. 
State-space equations S(t) = dx + Bu 
y(t) = Cy + Du 
Where A is called the state matrix, 
B is the input matrix 
Cis the output matrix 


D is the direct transmission matrix 


State-space representation 


1- If the forcing function does not involve derivative terms 


n n=l . 
Jy و‎ p alo +a, 11۲٢ لا‎ =u 
X, = X, 
let x, =y 
X, = X. 
X, y i 
, دد د‎ 


e 
ww 

| 
Sete 


X, Sut an X E — X, 
x وس چو‎ — EL. AL سوه سوه‎ — Ged, cw wu 
01 x = Ax + Bu 
yaly 

X, 0 1 0 0 | x, 0 
با‎ 0 0 1 0 | x, 0 a 
] ; IPM 0 X) 
; Jet ًِ l : "DP . 1+) 0 lu X4 
١ 0 y= 0 0... Of. 

0 0 0 1 0 

X, "Dun "ai "Ha? — à || X, 1 


Examples 


Consider the mechanical system shown in Figure `. We assume that the system is linear. The 
external force u(t) is the input to the system, and the displacement y(‘) of the mass is the output. 
The displacement y(z) is measured from the equilibrium position in the absence of the external 
force. This system is a single-input-single-output system. 

From the diagram, the system equation is 


my + by + ky =u 


This system is of second order. This means that the system involves two integrators. Let us define 
state variables x(t) and vc) as 


x(t) = y(t) 
l x(t) = SL) 
Then we obtain 


رند = پد 
1 1 
— ې ——(—KRy — bý)‏ = یا 
Sai ky y) -‏ 
=X "OR 3 4‏ ,× 
k b r : l‏ 
b uu‏ | 


‘The output equation is 
F 7344 


In a vector-matrix form, Equations ( 1 )and( 2 )can be written as 


1 0 i 4 0 
۳ = K b HM امه‎ 1 #١ 
X2 —— ر اا س‎ — 

m m nmi 





The sytem model can be written as 


wks‏ چ د 
otp"Tc FEO“‏ 


Then by defining state variables by 


X, = €g 
| X, = ê 
and the input and output variables by 
u = e, 
dan پس‎ 
we obtain 
0 1 
H = 1 R M TIG EN 
"| eet IER دل‎ 
and 


These two equations give a mathematical model of the system in state space. 


Example 2 





2- If the forcing function involves derivative terms 


y tay 


X, = y — Bou 
Xy =y- Bou - Bu =k, - Bu 


¥ = J > Pot - Bu ~ Bu = x, س‎ Bou 


(n—1) 
(n) VT. BE 17 4 d = ٨٨8 t bu T t A t d 


X = X 


* 


X, = 43 


(n) (n-1) ih 
8 +- mE ü 

e دص —- و‎ * 8 — (X 4 bau F bu n 

| (n1) — (0-2) ky, = "et 7 dr D 

(n1) ٣ n- 


y- pu - Bu =r- B,-2t B E DAN | X, 


X. = 0 1 l 0 ses 8 X1 Z 
' : : 1 T Q X5 
where Do, bi, fh», . .. , D, are determined from E د‎ | . i 
— l : 1 : + 
Po= be SH We = zs 
pim b, - abo 1 ae 0 0 0 1 Të CN 
rs | LM epa CE. “Hh. 404 Xp 
py = b, - ap, ¬ a, fio ١ Ion ۸ 
D E b. i dy» Hu à, p, 1 û4 Po | | 8 X2 
y-[1 0 0] + Bou 
x; 
p, — b, CARE th d, aD, i à, Bo 


~i 
CY 


Examples 


1-For the following system derive the state-system mode! 


Then ` M 


Es da ۸ b 


1 د‎ _ Kk b k DR 


and the output equation becomes 


OF 


and 


5 و ممصم 


Atzsslcse cari 


^d. 


Kb (NA KA d پا‎ 
a 


Hi |- 


it 


DEEN 


SEES 


EE 





wT ٢ دم‎ BS, 
dn idv er rr | : 
a". يم یت‎ i دوت‎ ai — Kiv — m} 
ea Let idt y 
a^ v ev , dn 
مسل ووو‎ oho + ky = bz c ku 
ar ade [HE 
h k 7 K 
E په پټ‎ E kee, Ach ks 
HH 1م‎ H 


with ihe standard 0 
jî + af + ayy = byt + bg Dali 


and identify a, az. Dn, bn, and b; as follows: 


D k - 0 I: 
e = "e KE جڪ‎ : diz, = 0 b, = DE H D- بو‎ 
dér Pray nmn déi 


2- Find the state-space for the following system u(s) T 





Y(s) | 160 (s + 4) 


^ wu q^ 2 
The corresponding transfer function is U(s) s TIN S T et s + 640 | 
Then, the corresponding differential equation is y -18 y - 192 y+ 640 y = 160 u + 640 u 





The state-space can be obtained as follows x, = y- Bou 
Let E = ym pü =p ü = $ = pu 

X, = y— يق اوس پل‎ Baue يل تار س ین‎ 
Where Bo B4 and DB, are determined from B. = 6, =0 


0 = b —a, f, = () 
p, =b, - a, fl, - a, B, =160 
f, = b,—a,f, — a, B, — a, f, =—2240 


Then the state equation is 


x. 0 1 0 E D, & 0 1 0 | x, 0 
ا‎ 0 0 L 1% #١) lu] ات ېل‎ 0 0 1 X i+; 160 D 
X5! Lë -a -ajx: |B; x,| |-640 -192 -18 دا‎ |—2240 
X X, 
y= |1 0 Of x, |+ fu DE 0 01: 


X3 e 


3 


Correlation between Transfer-function and 
State-space 

x(t) = Ax+ Bu 

y(t) = Cx + Du 


Let the state equation is given by 


Using Laplace transform technique yields 


sX(s)= AX(s)+ BU(s) 
(sI — A)X(s) = BUS) 


X(s)=(sI — A) BU(s) 

Y(s) =[C(sI — Ay B+ DIU(s) 
Ms) _ 
U(s) ` 
$(s) = (SI — A) ` and 


this 1s called State - transaition matrix 


G(s) = [C(sf - A)’ B +D]=Cġ(s)B + D 


Example :Find the Transfer function for the following system 


0 — [x 1 x 
i= , DN clu and y=(0 d' 
ve ` رر‎ Ol. ig 
L € 
Then ef — Al’ = = : 
و‎ +—S+— 
/ LC 
R 
od cht dl l 
S^ +—S+— 


Solving the time-invariant state equation 





1. Homogeneous sate equation 
: Scalar differential equation x = ax (1) 
The solution can be assumed as polynomial x(f) «5, FBSA br Ls. ED s. 
By substituting this into equation 1 yields 
b, 2b,t 13817 oe kot! +... = alb, + bt + bP +.........+ bf +...) (2) 
Equating the coefficients of the equal powers of t yields b, = ab, 
1 
b, = — = ah 
b, = —ab, = 1 a^b 
go LUE 342 0 
Lh 3 
kc x14 be 
Where b, can be determined by substituting t=0 into equation 2 
x(0) = b, 
Then, the solution is x(t) = ترو د‎ Puna tva +... kO 


=e“ xn 


e Vector matrix differential equation X = Ax 
Where x = n-vector A=n*n constant matrix 
If x is defined as 





2 | k 
X(t) = b, tat Ze Foses. Eb +... 
b, = Ab, 
bh, = Zb. = 1 43b. 
2 2 
b, = — Ab, = Ab, 
3*2 
D. = EH? 
k! 
m 1 77 1 Eak 
The solution is o =( 1+ 41+3 4ı "ایو‎ tt e pn 


-e^x(0)- d()x(0) ` 


Where ç(t) is called State-transition matrix 


Laplace transform approach to find the solution of 
homogeneous state equations 


۷ Ms sx(s) — x(0) = Ax(s) 


Taking the Laplace transform of both sides yields 
(sl — A)x(s) = x(0) 
Premultiplying both sides by ( SI — A) + yields X(s) = (sl - A) * x(0) 


2 
Where (sI — A) RE en 
S S A 


Then , take the inverse Laplace transform yields 
At At 
S [61-4 [=+ Att - ER T 


Therefore , the solution is 


x(t) = e^'x(0) = dat 


Properties of state-transition matrices 


For the time-invariant system 


g(0)=e4 =] 
b(t) =e" =(e“)' =p 
$^ (t) = d 
A(t, +t) = e t) = eth e* =A) = HA) 
P(t, —t)9(t, —1,) = Q(t —4) dt, —1)9 (t, — 1) 


(XO | = dnr) 


6- 


Examples 


1. Obtain the state-transition Q(t) and inverse state-transition matrix p-1(t) of the following 


x} | O 1 | x 
0 1 £||-2 3 
For this system A= 


mn و‎ =e” 2s fsr— An s-a- e sh i 
و‎ 3 
| s+ 3 s+1\s + DER fst 
IN ose GD دن‎ i [ene GC 2) ٣٢ 2) 
E چو‎ (3) 
Am y A 7 7 e 8 e? e uM aa 
-2e'42e" -e tie" 
2 igh d el 


x x 
Í ged Sg 
d H e d d'Be t —d 49 2t 


problem 


Find the solution of the control systems described by 


Y+5y+0y =J) 
where f(t)=0, y(0)=1, y(0)20 


y+1l0y+25y = f(t) 
where f(t)=0, y(O)=I1, y(0)=0 


x(t) = Ax + Bu 


2- non-Homogeneous sate equation 
y(t) = Cx + Du 


x(t) = e^ x(0) + | e^? Bu(r)dt = o(t)x(0) + | o(t—t)Bu(r)dt 
0 0 
Example : obtain the time response of the following system 
| To 1x] [o 
= T U 
Xs —— نه‎ ke 1 
Where u(t) is the unit step function applied at t = O. u(t) = 
--1 2t -1 2t 
م2‎ —e € —e 
At 
(t)-e = ره بر وهه‎ 
—2e د٢.6‎ -—e +2e 


STT ec Lie eg Lies سل‎ 
A , Jet? ے‎ eth) pL 7ے‎ 
[—T)= 
t $m )اب‎ ES ed 2\t- 
o [2e 2e e 27 
The response to the unit-step input is then D as 
| ۴ Zeite pene) gl) وج کي‎ 0 
x(f) sa e^'x(0) + j eh T) + je" He-r} -é g M77) A 2e 2M 77) IT 0: , 


OT 
- EY] [0.5 — PT + 0.5e 7 ] 
"S | 


“J TH i 
| ee 1 





In the special case where the initial state is zero, or x(0) = 0, the solution x(t) can be 
simplified to 


ھا = | 


Lxx(r) dëng 





|. ضحم‎ NHL SEES VMAS ES یا(‎ HATES هغ‎ mic د ماد تس‎ 


E 


Sylvester’s interpolation formula 
for compute ( e^*) 


The solution of el can be obtained by solving the following determinant equation 


LA à 2. AU ei 
1 A A» A, ev 


—] 7 
I EK a . ADU g” 


m mn 


e 4*5 s. A g^ 


By solving the above determinant e^* can be obtained in terms of the AK(k=0,1,2, m-1) 


Or it can be solved by writing the solution as 


e^ = a (f) a (D)4 a, (0)4 + ....+ a, (DAT 


Example 








e Compute e^t using Sylvester's interpolation formula if ol l | 
0 2 
LA g” 
lL A ead 
fy dx 
e — Where À are the eigenvalues of matrix A, and are the roots of the characteristic equation 
JAI — A| =0 
° Then 
Ao] fo ail, t 0 tj 
0 Al |8 wei | —2 2 Let 
I A چ‎ 


—— 
—— 





A = 
0 A+2 
A(A+2)=0 A=0, A,=-2 


() بل د اا د ar -2g = de‏ 


or ` Se >(4+21- Ae? 


mA 
لح‎ 
لاو‎ 
LG 
| 
E 
يم‎ 
ميټ‎ 
لن مسوم‎ ebe سی — دیل‎ 





An alternative approach is OL, (£) ta, DÉI zg 
a (f) + at (A, =e 
1 =0, 2 =2 
a,(t)=1 
a (t) - 2a, (t) =e" 


(0= 0-e) 


- e“ = a (ÐI ې‎ a, (0)4 = 1 + -e?)4 
1 21 
| صا‎ 
fi fem 
() سن‎ 


| D 1 
Problem compute e“ by two methods if A -| | 


Example 
Find the equilibrium points for the system described by 
y-ü-y)y-2y40.5y' =0 
Then evaluate the linearized Jacobian matrix at each equilibrium points. 


Determine the stability characteristics for the eigenvalues. 


Solution 
Find the state variable model 


let x, =y and x, = y gives the state variable model 


e = f(x) 


X5 2x, —0.5x° - (1+ x, )x, 


Equilibrium points are solutions of f(x) = 0, so each must have 


x, =0 and 2x, —0.5x; =0. 


0 2 —2 
The three solutions are — x, = 0 Xe =| gp and x,, = D 


The Jacobian matrix is 





0 1 

01 2 
W 2 

d. Ch 2 DEE DE so that 
of اح‎ its eigenvalues are +1 and -2 

OX 

3 
: [7 

a.p E: d eigenvalues are "ad j^ 

Ox |, -3 l 

01 1 dE 1 

—| = its eigenvalues are —+/,/!) 
E | Ü | : : ^ 


Matlab software 


To represent a polynomial in Matlab 
3 E 
p(s)=s +3s +4 


Define the polynomial as 

p=[ coefficients of the polynomial | 
p=[13 04]; 

To find the roots of the polynomial use 
R = roots( p) 


When the product of two polynomial is required the command ( CONV ) is used 
— Find the product of 


p=3s* +2s+1 
q=sS +4 


e p=[ 3 2 1],q=[1 4] 
e M=conv (p,q) 


DAP WN. 


How to represent the Transfer function in Matlab 


To represent the T.F in Matlab the following steps can be followed 
Define the numerator polynomial as 

Num=[ coefficients of numerator polynomial |; 

Define the denominator polynomial as 

Den=[coefficients of denominator polynomial ]; 

Sys = tf (num, den) 


10s^ +5s+1 


G(s) = ——— 
(s) و‎ +457 +25 +1 
Num=[10 5 1]; 
Den-[ 1 4 2 4]; 
Sys = tf(num,den) 
P=pole(sys) 


Z=zero(sys) 


Transformation from T.F to state-space or vice-versa 


257 +8546 
Uri ———— —— 
(s) x +85" +16 + 6 
-8 =]6 =6 | 
A4=| 1 0 0 B=|0 
0 1 0 0 
0-12 8 6| D=(|0| 
2 
+25+54 
ingen 
S د١٢5۵‎ 


The command to perform the conversion is 
[A,B,C,D]=tf2ss(num,den); 

Printsys(A,B,C,D) 

Or 

[num,den]=ss2tf(A,B,C,D,iu) 

% Where iu is specified for more than one input. 
Sys-tf(num,den) 

Let 


Numz[2 8 6]; den=[1 8 16 6]; 


For a system with multiple inputs.and multiple outputs; 


id 
Consider X, 0 L |x, ] Liz, 
= + 
Xx. —2) اه‎ | QO IL} u, 
»| |! O[x] [0 0 || ui 
y; 10 1 X5 0 Uu. 
A-[0 1; -25 -4]; 8-11 1; 0 1]; C=[1 0 ; 01]; D=[0 0;00]: 


[n1, d1]=ss2tf(A,B,C,D,1) ( to find transfer functions due to input 1) 
[n2,d2]=ss2tf(A,B,C,D,2) ( to find transfer functions due to input 2) 


(0 0 -25 
n2=0 1 5 d2=1 4 25 
Q 1 -25 
n11-[n1(1,1) n1(1,2) n1(1,3)] n11 =0 1 4 


n21-[n2(1,1) n2(1,2) n2(1,3)] n1220 0 -25 


n12-[n1(2,1) n1(2,2) n1(2,3)] 


n22-[n2(2,1) n2(2,2) n2(2,3) n21-01 5 
n22 = 0 1-25 


Sys11-tf(n11,d1,1) al ۹٩ 
U (s) s^-4s425 


Sys21=tf(n21,d2,1) hfs) . 25 


U (s) ei 45425 





Sys12=tf(n12,d1,2) Y (s) 1 $45 ^7 
U (ei s7 +4s+25 





Sys22=tf(n22,d2,2) bis) تعسو‎ 
U,(s) s «4s 25 


100 


Problems 


Find the state space model to levei-control system 





Find the Transfer function of the following system 


-3 1 0 
X + 1 
=a 0 1 


y-ll Oly 


Ta 


State-space representation forms 


— +b WY FDU 


The state space can be represented by four forms 
Let a system defined by 


Controllable canonical form 


1. 


This form is important in control system design. This form is 


0 0 X, 0 
1 0 x 0 
-+ u 
0 POI 0 
m ed 7 a X, 
Xi 
X) 
b, —a,b, | + byu 


0 
0 0 


X 


|o, WW a, D,‏ سو 


A state equation is 


Output state equation is 


d a, Dy 


— B. 


n-i 


e Observable canonical form 


This form is important for design an observer estimator . This has the following form 


A state equation is ^1 000..0 ca, " | b, 
is 1 0 0 0 --# ات‎ e P 
0 1 0 —É ر‎ 
= + 
- LU TEE RUE l E a 
1 (0 0 0 i =g .| x, 0 
. " 
Xa 
Output state equation is y= 0 0 1 
dod 
x 


UT SSSR ومو ومو‎ eres S 


3- Diagonal canonical form 


This can be derived only when the denominator polynomial involved only distinct real roots. 














Y(s) bys" څې‎ +....... +b لای‎ 
U(s) (s ٢ Xs + Dy sees (s+ nl 
C C C 
= by + 4 —— e... ې‎ 
d i S ۳ A S + P» Sor Py 
A state equation is 
x =p, 0 0 0 | x, 1 
X, 0 — D 0 0 X» 1 
E "ls پر‎ 
لا .-. .0 0 0 یه‎ Ix. 1 x, 
X, 0 0 E E CE " 
Output state equation is ysle e .... el «bw 


4- Jordan canonical form 


This form is used when there are multiple roots in the denominator. 

















let p, = رم‎ p, or Y(s) ` dys" FOS +... +5, 8 +5, 
UG) (s-p,J(s- p)... (s+ p, ) 
= by + i sda eh EE 
(s+p,) (stp) (stp) stp, s+ p, 
A state equation is irn l1 O0 ٥# 0 | x |0 
X; 0 -p O0 0 X» 
=p- A 0 | 
=| 0 U - 0 +} 1 o 
"NP 0 0 : Ó |x, 1 
D 0 0 0 -p, x, | 1 
a 
X) 
Output state equation is 
P 4 = lc, bs c, | + Dou 
X 
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Example 


Obtain state-space representations in the controllable canonical form, observable canonical 
form, and diagonal canonical form if the transfer function of the system is 


1 (ګ)‎ ` s+3 


U(s) WRC? or y+3y+2y=u+3u 


J+ ay + ayy = DU + Du + Dau 
n=2, a, =a,= 2,a,=3 b,-0, De 1 and رتا‎ 3 


Controllable canonical form 


x, 0 


سنا 
E.‏ 
© 


ERR 


x -0 -—a x ] 
2 2 1_2 ye) =B 17 


x1 
y(t) = D "och, b - a,b | x + Bou 


vf CES 
| UO 


شب شن E‏ نس 


2. observable canonical form 
a 0 -a | x, b, —a,b, Á | = 1 E d ^ | + 1 
= F u x | سک دت‎ 1 
i. 1 -a |x| |b -ab i 
; y(t) = [0 17 
y=|0 1]! |+bu ^2 

















X» 
3. Diagonal canonical form Y(s) M bat ha رن‎ urb, 
U (s) (s+ p, (s+ Pp») — (s-- p, ) 
=b, + و‎ S نې‎ ih 
SD SFP, روا‎ p, 
Hiel  s+3 و۸‎ +b 4+b, er? 
U(s) s?+3st+2 (s+p,\s+p,) (s+1)(s +2) 


Y) pa Go D. 2 , 
U(s) S--D, S+p, stl s+2 


Sete EM FG uh 


y - le RN y(t) = 2 4 











C2 
^d 


2 2 


= TE 


^ PDT 


problems 


1- Obtain the state-space representation for the following system in 


Y(s) ۹٥ 


— 


U(s) s +5s+6 


(a) controllable canonical form . 
(b) observable canonical form 


2- Obtain a state-space in a diagonal canonical form 


y+6y+lly+6y=6u | 


3- Find x1(t) and x2(t) of the system described by 


x} |0 1 |x, 
x,| |-3 IP 


Where the initial conditions are 


x(0)! | 1 
)د‎ |-1 
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Linearized using state variable 


Consider the following general nonlinear state variable model 


x= f(x,u,t) 
y=h(x,u,t) ۹ 
Assume a nominal solution X, (t) , u,(t) , and y,(t) is known . The difference between these 
nominal vector functions and some slightly functions x(t), u(t), and y(t) can be defined by 
ox = x(t) — x, (t) 
du = u(t) — u (t) 
óy = y(t) - y, (0) 


X, +ox = f (x, + Óx,u, + Ou,t) 


Then = f(x,,u,,t) + E4 Ox + pa ou + higher order term 
E u 


n 


y, toy = h(x, + óx,u, + Óu,t) 


= h(x,,u,,t)- 4 & + x ou + higher order term 
Ox |, Ou 


n 


Feature S EE‏ کسه :و سه هره 


: ELE 





Where [ ], means the derivative are evaluated on nominal solutions. Since the nominal 
solutions satisfy eq (1) , the first terms in Taylor series expansions cancel. 


For sufficiency small óx , ou and dy , the higher —order terms can be neglected, Then 


aed aZ 01 
Ox |, Ou |, 


EE a Z OU 
Ox |, Ou |, 


If x,(t) = x, = constant and if u,(t) = 0 =6,(t), then the stability of the equilibrium point x, is 


governed by 
Ox = E4 ÓX 
OF |, 


For this case , the system stability can be determine as following 


If all eigenvalues of the Jacobian matrix ( [6f/ 6x]) have negative real parts, the equilibrium 
point ( f(x)=0) is asymptotically stable for sufficiently small perturbation otherwise it will be 
unstable 


If one or more of the eigenvalues are on the jw axis and all others in the left-half plane, no 


conclusion about stability can be drawn from this model. 


Whether the actual behavior of the system is divergent or convergent will depend upon the 


neglected higher — order terms in the Taylor expansive . 


Except for the borderline jw axis case, stability of the nonlinear is the same as the 


linearized model at least near the equilibrium point . 
Equilibrium point is a point where the system can stay forever without moving 


Nonlinear systems frequently have more than one equilibrium point. 


Controllability 


Controllability isthe one of the important objective of the modern control. 


A system is said to be completely state controllable if 


For any initial time to , each initial state X( to) can be transferred to any 
final state X (tf) in a finite time tf > to by means of an unconstrained control UEL, 
Unconstrained control U(t) has no limit on the amplitudes. 


Consider the following system . It is clear that all initial state variable X, (t) can be controlled by the input 





signal u(t) 

d ag FU HX Ju -d3 "1 =o U à Ex, 

X, = Ux, کوھت‎ ke +u X4,|7| 0 — 5 Xs (E11 e 
X 0 0 en Ga D A 1 


A, = OX, FOR, —d x. tu 


KX, FEX, + Kx.‏ لر 


U POY encima ao ER on 


But for the system given below, the state variable X1 is not controlled by control u(t) . Then the system 


is uncontrollable , 





Ky X, ټی‎ D 0 || 0 
n x% |=| (0 -a 0 |x, |+} 1 lu 
T AN e x. 0 (NEE || x, 1 
i A | i x, 
i $n Ke y= |K, Koc Ax, 
s 5 


Example 2 


Determine whether the following system is controllable or not 
xj [-2 otsli 
ie +| lu 
X -1 —ljx, | |l 


| o2 
| 2 


controllability matrix is 


M, -|B 1= 


Rank of Mc =1 #n 
determinant =0 


Then the above system does not satisfied the controllability condition and it is 


uncontrollable. 


وو و ها سی ههجو صصاف etal‏ موق 


Method Il 


If the system is represented by canonical controllable form with non-repeated eigenvalues , then the 
controllability condition can be derived in the following 
Transfer the canonical controllable form into diagonal form by using the transformation 


X= Iz 

Then check the new matrix B". 
If matrix B" has no zero row, then the system is controllable 
This can be done by the following 

x= Ax+ Bu 

1z = ATz + Bu 

ż =T" ATz«T Bu = Az + B'u 

where A-T' AT, B'=T 8 


Where transformation matrix T can be found from 


DMS بر‎ 
دح‎ 
œ 


wë 771 JT 
2 S R e ^ 


4. 


Example use the transformation method state whether the following system is controllable or not. 


al [1 17x] fol : 
2 tou 
Xa =—2 l 3. 1 i 


The eigenvalues can be found by , ,À12jandA2--j 
Then AI i 4 e, 


4 ١ ۱ / 

A, A, J &J 
rel K 

Therefore (-)|ز2-‎ l 


!  rp-lp. 1 sc? 
B =T Bz-—— 
— AT) 1 


The B' has no zero row , then the system is completely state controllable 





Condition for complete state controllability in the s-plane 


A necessary and sufficient condition for complete state controllability is that no cancellation occur in 
the transfer function. 


Y) FFA 
Ex consider the following transfer function U(s) i (s 2.5\(s—1) 





Clearly, cancellation of factor ( s + 2.5) is occurred , then the system is uncontrollable. 


The state space form is * = 0 l ^ FS 7 
X. Zo cell xs 1 


controllability matrix is M. - [B اسا‎ i 
3 Ld 


Rank of Mc =1#n 
determinant =0 
Then the above system does not satisfied the controllability condition and it is 


uncontrollable. 


e m Rem m‏ شش مولام مش api‏ ته t hi‏ بج ومومو سهموذاب بو 


Output controllability 


In some cases, it is desired to control the output rather than the state of the system . 


A system is said to be output controllable if it is possible to construct an unconstrained vector U(t) that will 
transfer any given initial output y (t) to any final output y (t,) in a finite time ty St X tr. 
Any system described by state-space is completely output controllable if 


rank M,=|CA CAB CAB . . . C4" Di=m 


where m=no of rows of C matrix 


Where Mo is called output matrix 
Example : consider the system defined by 


The Mo is 
M,=|CB CAB Dļ=[0 1 0 


The rank of Mo 12 m 
Then the system is completely output controllable 
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How to obtain a controllable canonical form 
Steps to obtain the controllable canonical form 


The system must be completely state controllable. 


Find a characteristic equation 51 a A 1٢ EN -+ as" m +a bk ماب‎ a, 
n— 
Define a transformation matrix T as T= M, * Hi 
An- ed a 1 
a, a3 1 0 
W = 
a, 1 0 


منم 
Es‏ 
© 


Where ae are the coefficients of the characteristic equation 


|[ س 


A= 


SU 
| 


where 


لم 
Define a new state vector X‏ 


kal: 
x — Áx + Bu 
TX = ATX + Bu 


X -TATX +T” Bu = AX + Bu 


t 


EXAMPLE 


Consider the system defined by 


Transfer the system equation into controllable canonical form 
The controllability matrix is 
0 | -4 
M,=|B 4B 4’Bl=|0 o 1 
—3 9 
The rank of controllability matrix = 3 =n 
Thus the system is completely state controllable 


Characteristic equation is SI = Al — ei 4657 +1154+6 


00 6ک‎ ,a,_,=a,=l|l, a_,=a =6 


n—| 


ro Ar 


—- 
—— 


Then 


The inverse matrix is 


N بم‎ N 
tee 
|| 
EEN 
"i CO © 
O = CO 
= NÉI = 
اا ند د کل دس سن‎ 
por ( 
1-6 
= o T 
O © m 
EE 
|| 
~ 


po mc 
gel 0 C 


© بم D‏ 
on‏ تا =— 
[ ا کے 

|| 
EENEG, 
سم‎ O CO 
SÉ سم‎ C 
d Sg سم‎ 
a d 

| 

A) 

KN 

DY 

© 

= 


= MW 


The Transformation matrix T is 


Observability 


If the initial state vector , X(tO), can be determined from the observation of y(t) measured over 
a finite interval of time from tO [ tO x t < tf ]. the system is said to be completely observable. 


If all state variable has effect on the output , then the system is said to be observable. For 
example if the output equation for the diagonalized system with distinct eigenvalues is given by 


Y-Cx-[1 1 1]x 


Then the system is observable 





e If any state variable has no effect upon the output , then it can not evaluate this state- 
variable by observing the output. Then the system is called unobservable. 


e ifthe output equation is Y=[011]x then the system unobservable 





—G6 


e The concept of Observability is useful in solving the problem of reconstructing unmeasurable 
state variables from measurable variables in the minimum possible length of time. 


IND 
Cn 


RSS Seay بد هدر ېروم‎ Tan ate ai tear لم‎ attin ntm رس‎ e m دا دب‎ 


Condition for Observability 


The system is completely state observable if the rank of the following[ nmxn ] matrix is equal n. 


G 
CA 
The above matrix is called Observability matrix. M. = 
Or the determinant of the matrix z 0 0 
Example 1 
Consider the system given by can 
x O 1 O7% 0 
x15] 0. 0 FIS D 
X4 bo ose e x. 1 
d 0 5 1] e 
X3 


| -4 
State whether the above system is observable or Not 


From inspection ,it can be concluded that the system is unobservable . This conclusion is valid only for diagonalized 
system with distinct eigenvalues. i 


Observability matrix is C 
CA 0 5 1 
M,=| . |=|-4 -3 3 
-12 -13 -9 
CA" 


Rank of Mo =3 =n 
Then the above system satisfied the Observability condition and it is Observable. 





Example 2 
Determine if the system given below is observable 


u 





Solution 
The state and output equations for the system are 


The Observability matrix is 


C 3 4 
AM. an sums 
CA All ` — 16 


The determinant for this Observability matrix equals 0 . Thus, the Observability matrix does not have full 
rank, and the system is not observable 


E S 
P 
~i 


Condition for complete Observability in the ۶-06 


A necessary and sufficient condition for complete state Observability is that 
no cancellation occur in the transfer function. 


Consider the following system 


C(s) ` s+2 
R(s) i (s-- 2Xs 4-5) 


Show that if the cancellation is occurred the system is unobservable 


How to obtain a observable canonical form 


Steps 


The system must be completely state observable. 


Find a characteristic equati ES ess. SÉ 71-1 
quation کا‎ — A| =s” + as" ta, SHG, 
Define a transformation matrix Q as D= (w*M » 
a 0542 a, 1| 
HE o. 3E 2 x 2 د‎ 1٢ Û 
W = 
a, 0 0 
1 0 0 0 


Where a;s are the coefficients of the characteristic equation 


NJ 
dd 


4- Define a new state vector 


ei / 


x = Ox 

. H H Ab WM e = 
x= Áx + Bu (nn -à.. 
OX = AQX + Bu i 

1 E u where A= 
x = 0 AOX +O" Bu = Ax + Bu 

1 
0 I| -al 
b, —ü D, 
Dei MA 
B= C=CQ=|0 0 


b, — a,b, 


Example 


Consider a system defined by X = AX + Bu = l i ^ Ié 


-4 - 3 2 
y = CX = 1 Lx KE 
C 1 1 
Observability matrix is M, e CA = BECH 
Characteristic equation is SI — A| = 8 +28 +1 وح‎ + G8 + ds 


Then 


يي 


= — -— ee ZA 


The transformation matrix Q is 


ri E a 
Q = (W * M, Y A j|: jJ Ti 1 A-Q'AQ- E | 1 
Sigk | 


E V 
در‎ 
P 

"V 


In general the controllability and Observability can be illustrated graphically by 








Sco = completely controllable and completely Observable. 
Sc = Controllable but unobservable. 

So = Observable but uncontrollable. 

Su = Uncontrollable and unobservable 


problems 


For the following plants , determine the controllability . If the controllability can be determined by 
inspection, state that it can and then verify your conclusions using the controllability matrix 





d 


Determine whether the system given below is Observable or not. 


ál -2 ~i —-3]x1l [2 
x, |=| 0 -2 1 lx |+|1 llu] 
si عا‎ -—8.—9 x 12 
X 
y=|4 6 8] x, 
X3 


Determine whether or not each of the following systems is observable. 





Given the plant below , what relationship must be exist between C1 and C2 in order for the system to be 
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C? 





unobservable. 


2. 


r 


13 


5- Transfer the following system into the observable canonical form 





solved by analytical approach 


ily 


ich can not eas 


ind the response of complex systems wh 


Simulation 


is a method to solve and f 


Simulation 


Matlab simulink 


Matlab simulation ( software). Simulink has become the most widely used 
software package in academia and industry for modeling and simulating dynamic 


systems. 


Simulink is a software package for modeling, simulating, and analyzing dynamic 
systems. It supports linear and nonlinear systems, modeled in continuous time, 


sampled time. 


For modeling, Simulink provides a graphical user interface (GUI) for building 
models as block diagrams, using click-and-drag mouse operations. With this 


interface, you can draw the models easily. 


This is a far cry from previous simulation packages that require to formulate 


differential equations and difference equations in a language or program. 


mor hinotico. 


Simulink includes a comprehensive block library of sinks, sources, linear and nonlinear 


components, and connectors. 


. After a model is defined , it can simulate it, using a choice of integration methods, either 
from the Simulink menus or by entering commands in the MATLAB Command Window. 

The menus are particularly convenient for interactive work, while the command-line 
approach is very useful for running a batch of simulations. 


Using scopes and other display blocks, the simulation results can be seen while the 


simulation is running. 


In addition, the parameters can be changed and immediately see what happens. The 


simulation results can be put in the MATLAB workspace for post-processing and visualization. 


Examples 


How to model a simple system using simulink ? 

The model is to integrate a sine wave and displays the result along with the sine wave. 
To create the model, first enter simulink in the MATLAB command window. The Simulink Library Browser will 
appears. 

select the New Model button on the Library Browser's toolbar. 

Simulink opens a new model window. 

copy blocks into the model from the Simulink block libraries: 

The Sine Wave block from Sources library 

The Scope block from Sinks library. 

The Integrator block from Continuous library. 

The Mux block from Signals routing library. 


Then the system model is E DEI SE i 1 


integration of sine wave SCH 
pee nit: 


y" ann ږا‎ 


د مین مین رز اد 


وه 
Sine wav‏ 






Sine Wave 
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Integrator 


Effect of amplifier saturation 
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Example 2 


For the following system draw the system simulink diagram and find the system response for 


4 gli 
E | 


different values of system parameters. 


Bank angle, œ 






Plane of symmetry 


Bank angie, œ 
AU. "Ss 





: ری )ې‎ 
Pals} = O | 
Desired Bank angle 


bank angle 


2.Show the effects of disturbance on the system response. Assume step input 











Load 
E disturbance "TS 
Reference ` ` 2 | dÉ Ns) Actual 
speed | speed 
Rs) ioo? - 10 S I. L ei ` 
s? + 1405 + 100? . ئ).,ل‎ +1 | C S - +1 ھ‎ 
F lyweights Hydraulic . . Engine 


servo — | ۰ i 


3. Find the response of the following system under step input. 


Y(s) 








System identification 


System identification is the field of modeling dynamic systems from experimental data. 

in many cases the data may be corrupted by measurement noise and/or disturbances, which 
aggravates the modeling process. 

the techniques of system identification have a wide application area, 

The system is driven by input variables and disturbances .The output signals are variables 
which provide useful information about the system.. 


dynamic system can be described as in figure (1.1). 


Cistunmübaiuce 
wit} 


(Output 
yO 


Logout 
Wit} 





A dynamic system‏ ( 1 1ک لا 


EE 


Lay 


The basic information for System Identification 
Procedure 


The construction of a model from data involves 


three basic entities: 


— The input-output data 
— The set of candidate models 
— Determining the " best " model (model validation) 

















Parameter Medel 


validation | 






estimalic 





۸ تفسلقع‎ 
| kuowles]se | 





* 1. The input/ output data: 


The input-output data are sometimes recorded during a specifically 
designed identification experiment, where the user may determine which 
gp to measure and when to measure them and may also choose the 
input signals. The object of experiment design is thus to make these 
choices so that the data become maximally informative, subject to 
constraints that may be at hand. 

in other cases the user may not have the possibility to affect the 


experiment, but must use Data from Normal operation of the system 


6 


2. The set of candidate models: 


e a set of candidate models is obtained by specifying within which 
collection of models that are going to look for a suitable one. ` 

e This is the most important, and at the same time, the most difficult 
choice of the system identification procedure . 

e In many cases linear models are employed, without reference to the 
physical background of the system. 

e Such a models set, whose parameters used to adjust and fit the data 
and do not reflect physical considerations in the system, is called a 
black box model. 


3. Determining the "best" model in the set (Model validation). 


The assessment of model quality is typically based on how the models 
perform when they attempt to reproduce the measured data. 

The system identification procedure is a recursive one. It is quite likely that 
the model first obtained will not pass the model validation tests. then must 
be go back and revise the various steps of the procedure. 





Procedure of System Identification 


` € Design of Experiments 





|dentification of Model Structure‏ ظ 
Parameter Estimation‏ > 
Model Validation‏ * 


In the following figure the System Identification loop is shown: 


k | 
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The cycle can be itemized as follows 


Design an experiment and collect input-output data from the process to be identified. 


Examine the data. Polish it so as to remove trends and outliers, and select useful portions 


of the original data. Possibly apply filtering to enhance important frequency ranges. 


Select and define a model structure (a set of candidate system descriptions) within which a 


model is to be found. 


Compute the model parameters and then estimate the best model from the model structure 


according to the input-output data and a given criterion of fit. 
Examine the obtained model's properties by a validation data. 


If the model is good enough, then stop; otherwise go back to Step 3 to try another model 
set Possibly also try other estimation methods (Step 4) or work further on the input-output 


data (Steps 1 and 2). 
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Model Types 


Different approaches to system identification depending on model class 
Linear /Non-linear | E | 
a Parametric/Non-parametric 
e Non-parametric methods try to estimate a generic model (step responses, 
impulse responses, frequency responses). 
e Parametric methods estimate parameters in a user specified model (transfer 
functions, state-space method ). 


Model descriptions: 


e Transfer functions 
e State-space models 


Model Representation by Transfer-Functions: 


The most immediate way of parameterize G and H is to represent them as rational 
functions and let the parameters be the numerator and denominator coefficients. Such 
model structures are also known as black-box models 


*ARX-IMODEL Structure: 
(Auto Regressive with exogenous 
input) 


The simplest input-output relationship is obtained by describing it as a 
linear difference equation: 
y(k) +a y(k -1)*... a, y(k-n,) 
=bu(k-1)+...+b, u(k -n,) * e(k) 
Where: 
y(k -i)= Output in previous instant 
u(k -i)= Input in previous instant. 
n, = Number of a coefficients. 
n,- Number of b coefficients. 
e(k) = 1s the white-noise 


The white-noise term e(k) enters as a direct error in the difference 


equation. 


Let: 
O=[a,a,..a, bi 5, ... b, F 


A(q)=1+ aq” +... + gd" 


And 
B(q)=b,q" +...+ bog. 


Then G(,0) = | EE 


The above model is called the ARX-MODEL, where AR refers to the 
autoregressive part A(q)y(k) and X to the extra input B(q)u(X). 
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o ARMAX-MODEL Structure: 





(Auto Regressive Moving Average with exogenous input) 





"ch 





:OE-MODEL Structure: 
(Output Error model structure) 


bia) 
St u(k) + e(k) 





y(&)— 


e(k) 


. u(k) 





De (2.4) The OE model structure 


°BJ-MODEL Structure: 
(Box-Jenkins model structure) 


bly] Lig) 
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Finally the general model structure could be displayed as in following 
equation: 


| B(q) C(q) 
4 k)= k k 
(q)y(Kk) Fi; j+ ae ) 








And depending on which of five polynomials A, B, C, D, and F are used 
it may give rise to 32 different model sets. And the simpler models which 
already discussed are a subset of the general model structure, and 
possible by setting one or more of  A4(q),B(q),C(q),D(q),or F(q) 


polynomials equal to 1. 


خو رر راپ 


mane‏ یرم سه 
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Nonparametric Methods 


This section describes some nonparametric methods for system identification; the 
nonparametric methods could be seen as preliminary experiments. 


— Transient response analysis: 
= Step Response 


= [mpulse Response 


If the system which described by 


` Y(t)» G(q)u(t)+v(t) eg 1 
And 
U(t)-a 0 


zg t«O 


Is applied to Eq. (1) gives the output. 


y(t)-a Xg (k) +o | | (2) 


From this, estimates of g,(&) could be obtained as: 


a(t 240 Y€ D 3) 
0 


Which has an error ا(-٢-((‎ 


سا 
The transient performance of a control system is usually characterized by the use‏ 
of a unit step input.‏ 


And the typical performance criteria that are used to characterize the transient response to 


a unit step input are overshoot, delay time, rise time, settling time, static gain and time 
constant. 


٥. 


For example: 


For step response of the system shown in fig (3.1) we will calculate all the step response 
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fig(3.1) step response ٨ 
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This system clearly is a second order or higher under damped system because 
the system oscillates. 

The damping ratio of this system can be found using the percent overshoot. The 
percent overshoot can be measured of the step response by finding the ratio of 
the maximum achieved peak and the steady state output of the plot . 

And by using the matlab we found that: 


System: Sys 
| Peak amplitude: 3.72 | 
| Overshoot C3) 15.6 | بد‎ ep هم عمویعتها‎ 





| At time (Secr 3.59 unm NEN ! 
م -- سسرب‎ t 1 
ius LN m. | T— E 
AVE مها ماس‎ ie : ML Sys | 
e Y | Time (sec) 10.9 | pus : 
ووی‎ / ----------.4-.|Ampitude: té د‎ System: sys ; 
/ : (eee |! Ein : 
a - E - > : mal falc: Zt» c 
EE 
2 ub nix - می ابس‎ gem sy ee پر‎ eec cJ ee 
7 SEENEN (| System sys | 
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= i Time (secy: 14.8 
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fig(3.2) step response parameters Le? 


eFrequency response analysis: 


The input is a sinusoid, for a linear system in steady state the output will also be 


sinusoidal. The change in amplitude and phase will give the frequency response 
for the used frequency. 


For a discussion of frequency response analysis itis convenient to use the 
continuous time model. 


Y(s)= G(s) U(s) 


(4) 
If the input signal is a sinusoid 
u(t) = a sin(ax) (5) 
And the system is asymptotically stable, then in the steady state 
the output will become: 
y(t) = 0510) 604 +y) (6) 
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Where 


b=a|G(ja)| 
= arg| G( jo)| 


For example: 
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For the bode plot fig(3.3) we will calculate the frequency response parameters 
as follows: 


Estimating the Order 


The phase plot of the bode plot dips asymptotically to -180 degrees relative 
to the input. -180 degrees is two multiples of -90 degrees, so the system is 
at least second order. 

DC Gain 

The DC gain of a system can be calculated from the magnitude of the bode 
plot when s=0. 

The DC gain is: 

DC Gain = 10^-(M0/20) 

Where MO is the magnitude of the bode plot when jw=0. 

From the bode plot Fig(3.4): 

MO = 8.06 

DC Gain = 4 

The DC Gain of the system is 0.4 
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Natural Frequency ! 


The natural frequency of a second order system is the frequency when the 
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i System SYS 


Frequency (rad/secy 0.102 
Magnitude (dB): 8.06 


fig(3.4)Bode plot 
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phase of the response Is -90 degrees. 
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Experiments and Data Collection 


Often good to use a two-stage approach 
9 Preliminary experiments 


e Step/impulse response tests to get basic understanding 
01 system dynamics 


e Linearity, static gains, time delays, time constants, 
sampling interval 
® Data collection for model estimation 
e Carefully design experiment to enable good model fit 


e Operating point, input signal type, number of data - 
points to collect 


Preliminary Experiments — Step Response 


Useful for obtaining a qualitative information 
about the system 
۵ Dead time (delay) 
e Static gain 
e Time constants (rise time) 
e Resonance frequency 


Sample time can be determined from time constants 


(62 


Choosing the Right Structure 


Starting with non-parametric estimates (correlation analysis, 
spectral estimation ) 


® Gives information about model order and important frequency regions 


Prefilter input-output data to emphasize important frequency 
ranges 


^ Begin with ARX models 


Select model orders via 
e Cross-validation (simulate model and compare with new data) 
® Akiak's Information Criterion (AIC), i.e. pick the model that minimizes 


(1+ SE) 
\ NN ZC, 


t=1 
Where d is the number of estimated parameters in the model 
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Parametric Estimation Methods 


® Non-recursive/Batch (off-line) methods 
® Linear regression and (blockwise) least squares 
® Prediction error methods 
® Instrumental variable methods 
® Subspace methods (if possible} 

® Recursive (on-line) methods 
® Recursive Least Squares (RLS) methods 
۵ Recursive prediction error methods 
® Recursive instrumental variable methods 

. * Forgetting factor techniques and time-varying systems 

® identification methods 
® Blockwise sliding window least squares methods 


Principle of Model Validation 


® Compare model simulation/prediction with 
real data in time domain 

ə Compare estimated model's frequency 
response and spectral analysis result in 

frequency domain 

9 Perform statistical tests on prediction errors 


Keep Estimation and Validation Separate 


— Split data into two parts 
٥ one for estimation 
٥ one for validation 
—Apply input signal in validation data set to 
estimated model  : 
— Compare simulated output with output stored in 
validation data set کب‎ 
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Some commonly used input signals 


eThe input signal used in an identification 
experiment can have a significant influence on 
the resulting parameter estimation 
® Step function. 
® Pseudo Random Binary Sequence (PRBS). 
® With noise. 


9 Sum of sinusoidal. 
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PRBS. 


A pseudorandom binary sequence 
(PRBS) is a signal that shifts 
between two levels in a certain 
fashion. It can be generated by 
using shift registers for realizing a 
finite state system, and is a 
periodic signal. 


In most cases the period is chosen 
to be of the same order as the 
number of samples in the 
experiment, or larger 


When applying a PRBS, the user 
must select the two levels, the 
period and the clock period. The 
clock period is the minimal number 
of sampling interval 

Let u(t) be a PRBS that shifts 
between the values a and -a, and 
let it’s period be M, the 
autocorrelation function of PRBS is 
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White noise 


A sequence u(k) with 
autocorrelation function, 
zero except at rz à , 
having constant power 
spectral density function is 
known as white noise . 


The autocorrelation 
function of with noise is an 
impulse function 


Correlation Analysis 


The correlation analysis method that estimates 
the impulse response is useful only when the 
input signal u(k) is 0 zero-mean white 6 
signal. However, the input signal is not white 
noise in most real-world applications. The 
input u(k) and output y(k) signals therefore 
must be preconditioned before you apply the 
correlation analysis method. 
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Correlation Analysis 


The correlation analysis method uses the cross correlation between the 
input and output signals as an estimation of the impulse response, as 
shown by the followingsequation: 


> ulk —n)h(n) + elk) 





The input signal must be zero-mean white noise with a spectral density 
that is equally distributed across the whole frequency range. 


* Assuming the input u(k) of the system is a stationary, stochastic process 


and statistically independent of the disturbance e(k), the following 
equation is true. 


GE H (ON 


=8 


The above formula is called wiener-hopf equation. 


Correlation 5 


۰ Rat) represents the cross correlation function between the 
stimulus signal u(k) and the response signal y(k), as defined by the 


following equation. 
mepi 


) ylk +7) ulk) 
k-ü 
o Ru) represents the autocorrelation of the stimulus signal u(k), as 
defined by the following equation. ! N-:-1 
Rr) = : > u(k-t Tuk) 
IN La 
| k=0 
® Nis the number of data points. If the stimulus signal is a zero-mean 
white noise signal, the autocorrelation function reduces to the 


following equation. 
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Correlation Analysis 


where ouis the standard deviation of the stimulus white noise and 0(T) is the Dirac 
function. Substituting Ru(T) into the cross correlation function between the 
stimulus signal u(&) and the response signal y(k) yields the following equation. 


32.) B(k— D hlk) = o2h(2) 





There fore the impulse response h(k) is proportional to the cross correlation 
function between u(k) and y(k) 
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Correlation Analysis 


e If u(K) is PRBS signal. — Rt) =a tlt) 


® Where f is the clock pulse interval of PRBS. 
Using wiener-hopf equation 


R, (0) =a? ( Elk -0 hlk) = aho) 
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o Therefore the cross correlation function between u(k) and y(k) is 
proportional to the impulse response of the system h(k) 


Remark 


9 Correlation analysis has a main advantages over 
classical methods of step and sine wave testing 
methods in that the plant is kept at normal production. 

? The test input is superimposed with the operation input 

* The impulse response is unaffected by normal 
operating input and the noise, provided that these 


signal are uncorrelated from the test input 


Applications of the Impulse Response 


e The impulse response not only indicates the stability and causality of the 


system if feedback exists in the system, but also provides information on 


properties such as the damping, dominating time constant, and time delay. 


> Some of this information, such as the time delay, is useful for parametric 


model estimation. Therefore, you can use nonparametric impulse response 
estimation before parametric model estimation to help estimate the 


parameters. 


Spectral Analysis Method 


The frequency response analysis method of obtaining the frequency response Is 
straightforward but takes a long time to complete and is sensitive to noise. For 
these reasons, the spectral analysis method is used to estimate the frequency 
response function. 


You can use the spectral analysis method with any input signal. However, the 
frequency bandwidth of the input signal must cover the range of interest. 


Because the frequency response is the Fourier transform of the impulse response, 
applying the Fourier transform to both sides of the cross correlation function yields 
the following equation. 
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Spectral Analysis Method 


٥ G(e")is the frequency response of the system. 
e Oe") is the auto-spectral density of the stimulus signal ٨ 


e ,,(2) is the cross-spectral density between the stimulus 
signal u(k) and the response signal y(k). 


® You then can use the following equation to compute the 
frequency response 
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Applications of the Frequency Response 


e The frequency response gives the characteristics of the system in the frequency 
domain. You can use the frequency response to obtain useful information before 
applying parametric estimation. For example, you can use the frequency response 
to determine whether you need to pre-filter the signals or what the model order of 


the system is. 


eed 


* You also can use nonparametric frequency response to verify parametric model 
estimation results by comparing the frequency response of the parametric model 


with the nonparametric frequency response. 
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Examples 





SYSTEM IDENTIFICATION ofa first order system 


Step Response 


Consider the truck with a different load in 
the trailer, the driver experiences a slower 
acceleration as seen in the following step 
response due to increased mass and 
increased friction due to the added load 
bearing down on the tires. Assuming the 
truck has the same engine which produces 
20,000N of force 
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SYSTEM IDENTIFICATION ofa first order system 


This system appears to be a first order system, because the 


response does not oscillate and has a non-zero slope when t = 
0. For this reason we will model this system as a first order system. 


* DC GAIN 


The DC gain is the ratio of the steady state step response to the 
magnitude of a step input. The steady state step response can be 
determined from the plot of the step response, and the magnitude 
of the step input is given as u220000. 


SYSTE M IDENTI FICATION of a first order 


‘System: untitled 
| Time (sec): 348 | 
Amplitude: 33.2 


at 
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Step Response 


150 20 


Time (sec) 


100 





Am plitude 


system 


The DC Gain is: 

DC Gain = steady state output / step magnitude 
u = 20000; 

SS = 33.2; 

K = ss/u 

ka 
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SYSTEM IDENTI FICATION of a first order 





system 
SETTLING TIME 
To find the settling time of the system, 1 Step Response 
locate the time on the plot when the شت تس‎ 
l l e System. untitled! | 
magnitude crosses the desired percentage 30 EC N E 
1 a ^ System: untitled , Ampliude: 33:2 | 
of the final value. For instance to find the رر‎ Tine (sec). 132 
: : : 25 d ‘Amplitude: 29,9 | k 
10% settling time of this system, look for 8 d M | 
where the response reaches 0.9 of the / 
V i) 2 / 
final value: E / 
O20 C -» / 
ans = | 
H | 
29.8800 | | 
| 
ټمر اا‎ EE dd وح‎ M اچ‎ 
i H SQ 100 140 ZU لات‎ 3n] ann 
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SYSTEM IDENTIFICATION of a first order 


system 


9 TIME CONSTANT 


The time constant, t, of the system is the time at which the 


D 


response is 1 - 1/e = 63% of the final value. The relationship 
between the time constant and the pole of a system is: 


pole =-1/t. Other handy approximate relations for finding 
the time constant are: 





SYSTEM IDENTIFICATION of a first order 


system 


9 Because the 10% settling time is 132 seconds, 
the time constant is: t = 10% Settling Time / 
2.3 = 132s / 2.3 = 57s 


9 The time constant of this system is 
approximately 57 seconds. 


9 |t is also possible to find the time constant of 
the system by following the procedure for 
finding the settling time except for a 
percentage of 3796. 


SYSTE MI | D E NTI FICATI ON of the second order 


system 


UNDERDAMPED EXAMPLE 

Identify the system with the following Step Response 
response to a step input of magnitude 3. ` | 
ESTIMATING ORDER 
This system clearly is a second order or higher 
underdamped system because the system 
oscillates. The step response also has a zero 
initial slope when t-0, reinforcing the 
conclusion that the system is second order or 
higher. For this reason we will model this 
system as an underdamped second order 
System. 


Amplitude 





Time (sec) 
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SYSTEM IDENTIFICATION of the second order 


system 
DAMPING RATIO 


The damping ratio of this system can be found using the 

percent overshoot. The percent overshoot can be measured 

off the step response by finding the ratio of the maximum 
achieved peak and the steady state output Off the plot, the 14 


Step Response 
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/ System: 1 | 
Time (sec): 382 | 
p ‘System: uniledt | 
\ . Time (sec): 25 | 

i . Amplitude: 1 





Uh 


Time (sec) 


Amplitude 


percent overshoot is: 

“OS = 100% ( peak value - steady state output ) / steady 
state output 

peak = 1.39; 

SS = 1; 

os = 100*(peak-ss)/ss 

os= 39 

The damping ratio is: 

z = -In(%OS/100) / sqrt(pi?-1n?(960S/100)) 

dampingratio = -log(os/100)/sqrt(pi^2--(log(os/100))^2) 

dampingratio = 0.2871 

The damping ratio is 0.29, which is less than one, so this 
system is underdamped. 
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SYSTEM IDENTIFICATION of the second order 


system 


DC GAIN 
The DC gain is the ratio of the steady state step 


response to the magnitude of a step input. The 
Step Response 





steady state step response can be determined from 14 
the plot of the step response, and the magnitude of | f N 
the step input is given as u=3. The DC Gain is: 12 | \ (System: untied | 
DC Gain = steady state output / step magnitude | | | | Daten H 
ee? i \ e der 
ss = 1; | | Ke 
dcgain = ss/u , 8 
dcgain = E | 
$ DB | 
0.3333 | 
04 | 


The DC Gain of this system is 0.33. 
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system 


NATURAL FREQUENCY 

The natural frequency of an underdamped second order 
system can be found from the damped natural frequency 
which can be measured off the plot of the step response and 
the damping ratio which was calculated above. The natural 


Step Response 


frequency I5: System: untitled - 
Wn WA / sq rt(1 - z?) " | Time (sec): 3.82 | System: untiledt | 
| Time (sec) 113 | 
و‎ 108 — 


Amplitude: 1.39 | 
The damped natural frequency is: a oves 


Wy = 2p/Dt 

where Dt is the time interval between two consecutive 
peaks on the plot of the step response. 

dt = (11.3 - 3.82); 


wd = 2*pi/dt; 2 
wn = wd/sqrt(1-dampingratio^2) E 
wn = 0.8781 - 


The natural frequency of this system is approximately 0.88 
rad/s which means that in the absence of all damping, this 
system will oscillate at 0.88 rad/s. 
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order. 


degrees, so the system is at least second 


asymptotically to -180 degrees relative to the 


ESTIMATING ORDER 
The phase plot of the bode plot dips 


Plot. 
input. -180 degrees is two multiples of -90 
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Identify the following system with the Bod 


Bode Diagram 


system 


SYSTEM | DE NII FICATION of the second order 


SYSTE IM | DENTI FICATI ON of the second order 


system 


DC GAIN 

The DC Gain of a system can be calculated Bode Diagram 
from the magnitude of the bode plot when 
s=Q. 

The DC Gain is: 

DC Gain = 10M(0)/20 

where M(0) is the magnitude of the bode 
plot when jw=0. 

Using MATLAB: 

MO = -9.44: 

dcgain = 10^(M0/20) 

dcgain = 0.3373 

The DC Gain of the system is 0.34. 
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Frequency (rad/sec) 
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(bap) ased 


phase plot is at -90 degrees. 


the 


where w 56. is the frequency at which 
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(Bp) ه‎ ۲1251 


W on: 


response is -90 degrees. 
The natural frequency is: 
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The natural frequency of a second order 
system is the frequency when the phase of the 


Bode Diagram 


Natural Frequency 


system 
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SYSTEM IDENTIFICATION orthe second order 
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Bode Diagram 
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>) M IDENTIFICATION of the second order 


dcgain/(2*10^(M90/20)) 





The damping ratio of a system can be found with the 
DC Gain and the magnitude of the bode plot when 
the phase plot is -90 degrees. The damping ratio is: 
where M ss. is the magnitude of the bode plot whe 
The damping ratio of this system is 0.29 which is less 


than one, so this system is underdamped. 
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SYSTEM IDENTIFICATION of a high order 


system 
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Higher order systems tend to be difficult e 
to model because, in general, there are mems 7 
no analytical formulas to compute the | 
parameters of a model directly from the 
experimental data. However it is often 
possible to break down a modeling 
problem into the problems of deriving 
the parameters of first and second order 
systems. In this example, the Bode plot 
for the frequency response of a given 
set of experimental data will be 
analyzed similar to the previous 
examples but with unknown transfer 
function and steady state parameters 
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ECC DET 


En uU uU SES 
cata = [0.100200 1.0617 
0.101997 1.C0:147 - 1.17202:e-004H 
0.103912. — 1.C0:191 - 1.23347:e-024i 
0.105925 1.01236 - 1.302911e-024i 
0.107976 — 1.001283 - 1.30326:e-0041 
U.LLJUbS لال‎ 331 - 1. ¥34b5ve-UJ4L 
0.112202 1.001381 - 1.506223e-004i 
0.114376 ` 1.001434 - 1.581098e-004i 
0.116591 1.001468 - 1.659423e-004i 
0.118850 1.001544 - 1.741347e-004i 
0.121153 1,001602 - 1.827025e-004i 
0.123500 1.001663 - 1.916619e-004i 
0.125893 1.001725 - 2.010297e-004i 
0.128332 1.001790 - 2.108238e-004i 
0.130818 ^ 1.001858 - 2.210626e-004i 
0.133352 1.001928 - 2.317655e-004i 
0.135936 1.002001 - 2.429530e-004i 
0.138569 ^ 1.002076 - 2.546464e-004i 
0.141254 ` 1.002154 - 2.668680e-004i 
0.143990 1.002235 - 2.796413e-004i 






SYSTEM IDENTIFICATION of a high order system 


4945. 20254) 
GU, (IL — -U.LUUUZ7 - 1.52360Ce-UJ61 
AZSA. AJASE -0.101076 - 1.343327e-nl6i 
8413.8514l€ —-0.C00025 - 1.362736e-036i 
8576,95398€ — -0.C00024 - 1.28561^e-026i 
6743.124580 -0.000023 - I.Z214797e-0061 
6912.509361 -0.000023 - 1.146956e-006i 
9085.175757 -0.000022 - 1.082901e-006i 
9261.187281 -0.000021 - 1.022419e-006i 
9440.608763 -0.000020 - 9.653114e-007i 
9623.506264 -0.000019 - 9.113908e-0071i 
9809.947127 -0.000019 - 8.604793e-007i 
10000.000000 -0.000018 - 8,124093e-007i]; 


freq = data(:,1); 
resp = data(:,2); 


frd model = frd(resp, freq); 


figure (1); 





The data must now be turned and stored 
as a Frequency Response Data (FRD) 
model. To store the data as a FRD model, 
the frequency and response must be 
stored separately with the following 
commands” 

freq = data(:,1); 

resp = data(:,2); 

Now that the frequency and 
corresponding response are separate, 
they can be stored using the FRD 
command. 

frd model = frd(resp, freq); 


| SYSTEM IDENTIFICATION of a high order system 


With this model we will use the Bode plot of the 
FRD model to create and compare a second 
order system and its corresponding Bode plot. 
The following will be the resulting Bode plot. 
Notice on the above magnitude plot the fairly 
level graph from 0.1 to 3 followed by a peak and ` CG 8٢ Es 
adrop. At high frequencies, the magnitude pun - B 
graph decreases at about 40 dB/dec. Also, | 2 
notice that the phase plot decreases from O to - 
180. These are all characteristics of a second 
order system. So to start the modeling process, 
construct a new second order transfer function 
to capture these features of the Bode plot. 








SYSTEM IDENTIFICATION of a high order system 


To begin we will use the following second order 
transfer function model: 


G(s) in Kae 


Tus 


, BEER 
im یت‎ re m m re e e 


| +} Figure No. 1 
iiem 





1 
2 : 
=} +o) 44 
W Qj, 


Ti 
Using the graph we can locate the frequency, Wo) 
and magnitude, Mp, of the peak and the DC 
gain, kae By dragging the cursor along the plot 
of the graph, one can find the different values 
and frequencies. Find the magnitude at the 
lowest frequency of the graph which will be a 
good approximation of the DC gain. Also find 
the frequency and magnitude of the peak as 


shown. 


d Frequency (radésec): 2.54 
| Magnitude (dB): 2.45 


g‏ سو مس وم e at‏ مهه Re‏ خو وو 


System: exp frd model 


SYSTEM IDENTIFICATION of a high order system 


The relationship between features of the Bode plot 
and the parameters of the second order system are as 
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Where w, is the natural frequency where the graph 


p? 254 





i.e. 2.45 dB). 


begins to turn and C is the Gamping ratio. 
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From the graph d 


rad/sec, and Mp = 1.33 
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SYSTEM IDENTIFICATION of a high order system 


These differences suggest adding a so called lead SE 
compensator transfer function to the second order ie =e IU UE 1339 De 


model. A lead compensator has,.the following form: 
1 


The lead compensator has a zero at -z and a pole at - 
p, which are defined by: 





& = W 


als 


Here, is the amount of phase at the phase peak 


and is the frequency of the phase peak. ` 
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SYSTEM IDENTIFICATION of a high order system 


EES Ee‏ وغم دا 


The next step of the system identification 
process is to find the frequency at which the 
phase difference between the FRD model and 
the second order model is greatest. This 
frequency will be the of the lead 
compensator, and the&tifference in phase will 
bethe ofthe lead compensator. 
A good estimate for these values would be: 

= 30%ad/sec, = (-91) - (-175) = 84 degrees. 
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SYSTEM IDENTIFICATION of a high order system 


T———— EE وم وه‎ MÀ 


ure No. 1 





Then graph the Bode plot the of the second 
order transfer function multiplied with the 
lead compensator as shown. Notice that the 
peak is higher than expected. This is because 
of the influence of the lead 

compensator. Use just the resulting values of 
6 =0.4127 and ٧۷ 3.1282 without the 
equations as a inital values for further 


iterations. 


| Frequency irad/secy: 2.52 
| Magnitude (387 267 


Selen frd model | 
Frequency (rad/sec). 252 
| Magnitude (dB): 245 | 





A good starting point for a iteration would be to find 
the difference between peaks of the FRD model and 
higher order model and multiply that value to 

¢. Remember to convert from dB. Here the 
difference it found to be 6.11 dB (i.e. 

2.02). Multiplying zeta by this value gives the 
following plots. 

the peaks are nearly coincident and may be a 
sufficient model for many cases. Often times, 
matching the lower frequency values where the 
magnitude is highest is more important then 
matching the high ones where the magnitude drops 
off rapidly. Proceeding two iteration as follows ` 
results in an even better graph. 





SYSTEM IDENTIFICATION of a high order system 


Above, the peaks are nearly coincident and may be 
a sufficient model for many cases. Often times, 
matching the lower frequency values where the 
magnitude is highest is more important then 
matching the high ones where the magnitude 


drops off rapidly. Proceeding two iteration as 
follows results in an even better graph. 








SYSTEM IDENTIFICATION ofa high order system 


e % Iteration 1 


zeta = zeta * 10^((8.56-2.45)/20) 
second tf = 1/((s/wn)^2 + 2*zeta*(s/wn) + 1); 


% Iteration 2 

wlead = 35; 

philead = 81 * (pi / 180); 

a = (1 - sin(philead))/(1+sin(philead)); 
zlead = wlead * sqrt(a); 

plead = wlead / sart(a); 

leadcomp = (s/zlead + 1)/(s/plead + 1); 


96 ...Iteration X 

zeta = zeta / 10^((3.0)/20) | 
second tf = 1/((s/wn)^2 + 2*zeta*(s/wn) + 1); 
figure(1); 

bode(frd model, second tf * leadcomp,{1, 1011: 
title(highermodel’); 
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summary 


Nonparametric model estimation is simple and more efficient, but often less accurate, than 
parametric estimation. 

you can use a nonparametric model estimation method to obtain useful information about a 
system before applying parametric model estimation. For example, you can use 
nonparametric model estimation to determine whether the system requires preconditioning, 
what the time delay of the system is, what model order to select, and so on. 

You also can use nonparametric model estimation to verify parametric models. For example, 
you can compare the Bode plot of a parametric model with the frequency response of the 


nonparametric model 


E? 
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